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We report on development of a new feature of the code HFODD, allowing for the 
angular-momentum projection of cranked symmetry-unrestricted Slater determinants. 
After a brief overview of the main theoretical building blocks and formalism, we present 
several preliminary applications. In particular, we discuss the case of a well-deformed 
rotational band in 156 Gd, and we show the emergence of uncompensated poles in the 
overlap kernels calculated in an odd- A nucleus 155 Eu. 



1. Introduction 

Self-consistent mean-field method is practically the only formalism allowing for 
large-scale computations in heavy, open-shell nuclei with many valence particles. 
Inherent to the mean-field method is a mechanism of the spontaneous symmetry 
breaking, which is essentially the only way allowing to incorporate a significant part 
of many-body correlations into a single intrinsic (symmetry-violating) Slater deter- 
minant. Nuclear deformation and the emergence of rotational bands after applying 
simple cranked mean- field approach, i.e., after an approximate restoration of rota- 
tional invariance, is one of the most spectacular and most intuitive manifestations 
of the spontaneous symmetry breaking in nuclei. 

An increasing quality of spectroscopic information demands more and more ac- 
curate theory, pushing the efforts beyond the mean-field approximation towards, in 
particular, symmetry-conserving formalisms and symmetry restoration. 1 ' 2,3 Hereby, 



1 



9, 2008 0:47 WSPC/INSTRUCTION FILE rzut!56-12 



2 H. Zdunczuk, J. Dobaczewski, W. Satula 

we report on development of a new theoretical tool allowing for the angular- 
momentum projection (AMP) of cranked symmetry-unrestricted Slater determi- 
nants. 

The standard method used to develop symmetry-conserving theory, starting 
from intrinsic (deformed) wave function |<£>), obtained within the mean- field ap- 
proach, is provided by the projection techniques onto eigenspaces of symmetry 
operators. There are two practical realizations of the projection methods: more 
fundamental and elaborate variation after projection (VAP) and slightly less ad- 
vanced projection after variation (PAV). In the past, many calculations have been 
performed within the angular-momentum PAV method, where non-rotating states 
have been projected, see, e.g., Refs. 4 ' 5,6 and the reviews in Refs. 2 ' 3 . The crank- 
ing method provides the first-order approximation to the angular-momentum VAP 
method. 1 However, after the ground-breaking studies in Refs. 7 ' 8 , calculations based 
on the AMP of cranked states have not been performed. Here, we present the first 
results of our recently developed AMP method of cranked symmetry-unrestricted 
Hartree-Fock states. The procedure we use has been implemented within the code 
HFODD (v2.25b). 9 < 10 < n < 12 

We determine the optimal product wave function |$) by using the cranked self- 
consistent Skyrme-Hartree-Fock (SHF) method. The Ritz variational principle: 

(*\H-u,I y \*) 

($|$) u ' w 

is equivalent to a variation of a local energy density functional (LEDF) with a 
supplementary constraint on the average value of the projection of angular momen- 
tum on the axis perpendicular to the symmetry axis (y in our case). The variation 
of LEDF in the isoscalar-isovector t = 0, 1 representation (potential part) can be 
formally written as: 

§s (Skyrme) = g ^ f ^ [ n [ TE \ r ) + H f°\ r )], ( 2 ) 

t=o,i 



where 7i t (r) is expressed as a bilinear form of the time-even (TE) p, r, J and time- 
odd (TO) s, T, j local densities and currents, and by their derivatives, 13 

n[ T E) (r) = C p tP 2 t + C? p p t Ap t + ClfHT t + C(j 2 t + C7 J Pt V • J t , (3) 

Wf 0) (r) = C t s s t 2 + C t As s t As t + Cjs t ■ T t + C\f t + C7>s t ■ (V x j t ). (4) 

By taking an expectation value of the Skyrme force over the Slater determinant, 
one obtains the LEDF (3)-(4) with 20 coupling constants C that are expressed 
uniquely through the 10 parameters Xi,ti,i — 0,1,2,3, and W,a of the standard 
Skyrme force. Because of the local gauge invariance of the Skyrme force, only 14 
coupling constants C are independent quantities. The local gauge invariance links 
three pairs of time-even and time-odd constants in the following way: 

ci = -ci, c! = -cf, c?J = c7 J . (5) 
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2. Angular-Momentum Projection 

A deformed solution |$) of Eq. (1) is a superposition of eigenstates of the angular- 
momentum operator. The angular-momentum conserving wave function can be ob- 
tained from the state |$) by applying the angular- momentum projector: 

\IMK)=P I MK \<S>). (6) 

The AMP operator P\ AK is a projector onto angular momentum / with projection 
M along the laboratory z axis, and reads: 14,1 

P mk = J R{tt) dSl, (7) 

where f2 represents the set of three Euler angles (a, f3, 7), Dj^ K (Q) are the Wigncr 
functions, 15 and R(Q) = e ~ laIz e~ ll3I ve~ ljIz is the rotation operator. 

The state \IMK) is no longer a product wave function, but a complicated 
superposition of Slater determinants. The operator P[f K is not a projector in the 
mathematical sense. 1 It extracts from the intrinsic wave function the component 
with a projection K along the intrinsic z axis of the nucleus. Since K is not a good 
quantum number, all these components must be mixed, and the mixing coefficients 

(i) 

g K must be determined by the minimization of energy. This i^-mixing is taken into 
account by assuming the following form for the eigenstates: 

\IM)^ = Y j $\IMK)^Y.9k P mk\*>), (8) 

K K 

and by solving the following Hill- Wheeler (HW) equation: 

Hg = EMg, (9) 

where Hk'k = (®\HPk>k\®) an d J^K'K = (®\Pk'k\®) denote the Hamiltonian 
and overlap kernels, respectively. 

3. Application to 156 Gd 

We solve the cranked SHF equations for 156 Gd by using the code HFODD for the 
SIII Skyrme-force parameters 16 and spherical harmonic-oscillator basis composed 
of A^o=10 shells. Then, we calculate the Hamiltonian and overlap kernels by using 
50 Gauss-Chebyshev integration points in the a and 7 directions and 50 Gauss- 
Legendre points in the [3 direction. 12 Figure 1 shows probability distributions Wi 
of even angular- momentum components I in the intrinsic crankcd-SHF states \&i y ) 
constrained to (I y ) = I y . The probability of finding the I component can be calcu- 
lated from: 1 

Wl=^I v \Pk K \^Iy)- (10) 

K 

The curves correspond to cranking wave functions with averaged angular momenta 
(I y ) = 0,8,12,16, and 20h. One can see that for low angular momenta (e.g., for 



9, 2008 0:47 WSPC/INSTRUCTION FILE rzut!56-12 



re 

o 
a. 




0.00 



4 8 12 16 20 

Angular momentum I [ ft ] 



Fig. 1. Probability distributions of even angular-momentum components I in 156 Gd. 



I y = 0) the maxima of the distributions do not lie near the same values of I ~ (I y ). 
Similar results have been obtained by Islam et al. 7 and Bayc et al. s . 

In Fig. 2 we show similar probability distributions Wk of the projections K, 



Wk = (*i v \PIck\*i v )/Wi, 



(11) 



projected from the state with (I y ) = 10h in 156 Gd. One can see that the K = 
component dominates for all angular momenta, while the K > components 
increase with the increasing angular momentum. The K — 1 component is the 
second in magnitude after K = 0, which illustrates the build-up of the Coriolis 
coupling in a rotating intrinsic state. Note that only even values of K can be 
projected from the non-rotating (I y ) = state, while all K > appear in the 
cranked state. 

In order to find energies of the AMP states, we solve Eq. (9) by diagonalizing 
first the norm matrix: 



J\fr] = nr]. 



(12) 



The non-zero eigenvalues (n m ^ 0) of TV are used afterwards to built the so-called 
natural states 



\m) = -^=J2vln ) \IMK), 



(13) 



K 
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Fig. 2. Probability distributions of projections K in the even angular-momentum components / 
projected from the state with (I v ) = 10H in lo6 Gd. 

that span the subspace called collective subspace. Final diagonalization of the HW 
equation (9) is performed within the collective subspace. In this subspace the prob- 
lem reduces to the standard hermitian eigenvalue problem. 

In practical numerical applications we use the cut-off parameter £ to construct 
the collective subspace, by keeping only the norm eigenstates (12) with n m > (. The 
test depicted in Fig. 3 shows the stability of projected solutions with respect to the 
number of the norm eigenstates kept in the collective subspace. The analysis was 
conducted for angular-momentum < I < 20 states projected from the deformed 
solution \&i v= io) obtained by solving the cranked SHF equations with the constraint 
(®i v =io\Iy\®iy=w) = 10ft- The test clearly shows that, starting from a certain point, 
the obtained solutions are stable (plateau condition) . Only for very small values of 
£ < 10~ 8 the method becomes numerically unstable. 

Figure 4 shows rotational bands (b)-(d) calculated in 156 Gd in comparison with 
the experimental data (a). Band (b) represents the average mean- field energies 
obtained from the cranked SHF calculations by constraining solutions to (I y ) = I. 
In bands (c) and (d), we show energies obtained by the AMP from the (I y ) = 
and (I y ) = I, respectively. We see, that band (c) is much higher than bands (b) 
and (d), which indicates that the PAV from the (I y ) =0 state is not an adequate 
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Fig. 3. Dependence of energies of projected states on the number of the norm eigenstates kept in 
the collective subspacc. Angular momenta < / < 20 were projected from the state having the 
average projection of angular momentum (I y ) = 10h. 



method of describing nuclear rotation. Note that the mean-field energies (b) and 
the AMP energies (d) corresponding to (I y ) = I turn out to be very similar to one 
another. This shows that the AMP from the (I y ) = I states constitutes a correct 
symmetry restoration of the approximate VAP solutions realized by the cranking 
procedure. The remaining discrepancy with experimental data must be attributed 
to pairing correlations, which arc not included in our SHF solutions. 

When calculating the Hamiltonian kernels within the LEDF approach, one has 
to use transition density matrices between rotated states, 

(<t>\a+a a R(n)\f>) 

Pafiity = —~ , (14) 

which arc singular whenever the rotated states are orthogonal. In particular, the 
transition multipole moments, 

QXfiiP-) = ^(QA^Wa^), (15) 

a/3 

become singular for certain values of the Euler angles fi. This is illustrated in 
Fig. 5, which shows absolute values of the neutron and proton overlap kernels in 
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Fig. 4. Rotational bands in 156 Gd nucleus: the four bands represent: (a) experimental data, (b) 
cranked SHF calculations, (c) band projected from the state |<E>/ y= o), and (d) band projected 
from the state \& Iy=I ) into /, for I = 0, 2, . . . 20. 

155 Eu (left panel), and the transition quadrupole moments Q2o(^) in 155 Eu and 
156 Gd (right panel). Calculations have been performed for axial shapes of nuclei, 
for which only the rotation about the Euler angle (3 matters. The neutron overlap 
kernel corresponds to an even number of particles (N = 92), and is always positive, 
although at (3 — 90° it becomes as small as 10~ 14 . On the other hand, the proton 
overlap kernel corresponds to an odd number of particles (Z — 63), and it three 
times changes the sign in the interval of 0° < (3 < 180°. Consequently, the transition 
quadrupole moment of 156 Gd is a regular function, while that of 155 Eu has three 
poles. 

Of course, when calculating the matrix elements of multipole operators between 
the rotated states, the transition matrix elements (15) are multiplied by the over- 
lap kernels and the poles disappear. However, such a compensation is absent for 
kernels corresponding to higher powers of densities, viz. for the density-dependent 
terms of the Skyrme interactions, or for the direct-Coulomb-energy terms, or for 
the exchange- Coulomb-energy terms in the Slater approximation. How to treat such 
singular kernels within the AMP methods is currently an open and unsolved prob- 
lem, similarly as is the case for the particle-number-projection methods recently 
discussed in Ref. 17 . 

4. Conclusions 

We have developed a new theoretical tool allowing for the angular-momentum pro- 
jection of cranked symmetry-unrestricted Slater determinants. In the present work, 
we presented preliminary applications pertaining to rotational bands in a well- 
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Fig. 5. Left panel: Dependence of the absolute value of the overlap between rotated states on 
the Euler angle /3 in 155 Eu. Solid line shows the neutron overlap that is always positive. Dashed 
and dotted lines show these segments of the proton overlap where it is positive and negative, 
respectively. Right panel: Transition quadrupole moments Q20 versus (3 in 156 Gd (thick line) and 
155 Eu (thin line). 

deformed nucleus 156 Gd. We showed that energies projected from non-rotating and 
rotating states are significantly different and similar, respectively, when compared 
to the cranking mean-field energies. We also showed that, due to vanishing overlaps 
between rotated states, transition energy kernels may have uncompensated poles 
in function of the Euler angles. 
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